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Abstract 

Theoretical calculations show that above 80 km in the earth’s atmosphere the production of 
vibrationally excited ozone by chemical processes leads to number densities which are usually larger 
than those expected for local thermodynamic equilibrium. Quenching of hi^ly excited molecules 
produced in O-^O^-t-M -*■ 03't-M provides a significant source of the lower lying states above the 
mesopause while the 9.6m emission of 03(0,0, 1 ) is a miyor sink. Analysis of available laboratory 
results implies that reactions involving excited ozone play a significant role in the global ozone 
balance despite the relatively small abundance of the molecule. However, this effect is implicit in 
many of the rate coefficients currently used in stratospheric calculations. In the upper mesosphere 
and lower themosphere, where the excited state populations differ from those for thermal equilibrium, 
published reaction rate data are not necessarily applicable to aeronomic calculations. 
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THE AERONOMY OF VIB RATIONALLY EXCITED OZONE 


Introduction 

The agronomic importance of electronically excited species is well-known. Prominent exampkis 
are reactions of the metastable oxygen atom, 0(* DI+HjO -*> OH+OH and 0(^ D)+N 20 -*■ NO+NO, 
which provide stratospheric sources of odd hydrogen and odd nitrogen, while excited states of 
molecules such as 02 (^ Ag) and OH(v' > 0 ) lead to airglow emissions that act as remote sensors of 
conditions in the mesosphere and lower thermosphere (LleweUyn et al., 1973 ; Noxon, 1978 ; 
Frederick et aL, 1978 ; Frederick, 1979), Cvetanovic ( 1974 ) has given an overview of the species 
mentioned above and their role in upper atmospheric processes. Recently, attention has focused 
on the potential significance of electronically excited ozone even though there is doubt concerning 
the very existence of bound excited states below the dissociation limit. (Hay et aL, 1973 , 1975 ; 
Wraight, 1977 ; Prasad and Burton, 1979 ). On the other hand, experimental evidence for the chem- 
ical production of vibrationally excited ozone in the atmosphere appears definitive (von Rosenberg 
and Trainor, 1974 , 1975 ). This paper combines the available data pertaining to production and 
loss of vibrationally excited ozone so as to estimate its abundance in the earth’s atmosphere. Given 
this infonnation, together with oublished laboratory results concerning the influence of vibrational 
excitation on reaction rate coefficients (Gordon and Lin, 1973 ; Kurylo et aL, 1975 ), one can esti- 
mate the potential significance of excited O 3 in the chemical procei^s which occur in the upper 
atmosphere. 

Aeronomic Significance of Excited States 

In a number of known cases, reactions involving excited species proceed more rapidly than 
when the reactants are in the ground state. Specific examples are the processes: 

Of^Pl + Oj-Oj +O 2 kj (Rl) 

0(‘D) + 03 -^ O 2 +O 2 k 2 , (R2a) 

-O 2 t-0(3p) + 0 (^P) k 2 b (R 2 b) 
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where kj = 1.5 x 10"^^ exp(-2218/T) cm^s“^ but k2 * ^2,+ k2b ■ 2.4 x cm^s“* , and 

N(^S) + 02 "►NO + O2 kj (R3) 

N(2d) + O 2 no + Oj k4 (R4) 

with k3 ■ 4.4 X 10"^^ exp(-3220/T) and k4 ■ 6 x 10“*^. All rate coefficients listed above are as 
given by the NASA Panel (1979) except for k4 which is from Lin and Kaufman (1971). The extra 
1.97 eV available in R2 over that in R1 lends to a gas kinetic rate coefficient. In the case of R4, the 
activation energy of R3 is absent but the pre>exponential term is essentially unaltered. 

The efficiency of a reaction depends on details of molecular scale interactions and one can- 
not necessarily assume that additional vibrational and electronic energy enhances the rate coeffident. 
Nonetheless, the above reactions clearly show that there are cases where exdted states are aero- 
nomically relevant. A process of particular importance to the stratosphere is the reaction of nitric 
oxide with ozone. Studies of the reaction; 

NO + 03(0, 1 .0) NO2 + Oj kj (R5) 

by Gordon and Lin (1973) and Kuiylo et al. (1974) show k^ to be a factor of 22 larger than the 
rate coefficient, kg of: 

NO + 03(0,0,0) - NO2 + O2 kg (R6) 

at a temperature of 350° K. This is approximately the relationship one would expect of the vibra- 
tional energy of 03(0, 1 ,0) were used to overcome the activation energy of R6, that is kg kg 
exp [hcK010)/kT] where v is the term value of the exdted state. In the presence of a Boltzmann 
distribution of vibrational levels the exponential scaling of rate coeffidents, which in this case has 
empirical justification, then implies; 

kg [03(0,0,0)] [NO] - kg [03(0, 1,0)] [NO] 

The rate of odd oxygen loss in reactions of NO with ground state ozone is equal to that when the 
excited (0,1.0) state is involved. Hence, despite the small number densities, reactions involving ex- 
cited vibrational states of ozone can play an important role in certain upper atmospheric chemical 
processes. If the rate coefficients for the reaction of NO with 03(0,0.!) and OgC 1,0.0) scale as docs 


that for 03(0, 1,0), then RS is of minor importance compared to reactions involving excited ozone. 
Similarly, studies of R1 with 03(0,0, 1) as a reactant indicate a room temperatture rate coefficient 
ISO to ISOO times faster tl^an when ground state ozone is involved provided the process proceeds 
by reaction producing O2 rather than vibrational deactivation (West et al., 1976). 

Given that reactions of excited ozone are significant in some circumstances, it becomes 
important to consider that published rate coefficients refer to the particular vibrational distributions 
present in the laboratory. If one considers the series of processes: 

X + 03(0 XO + Oj kj 

where i refers to a specific state, the rate coefficient determined in the laboratory, k, is: 

ic » Zkj 

i ‘ 

Only a limited number of terms contribute to the sum since kj can increase with i to at most a gas 
kinetic value while [03(1)] / [O3 ] decreases to a negligible value for sufficiently large i. If the 
fractional population of the excited states are identical in the laboratory and atmosphere, then the 
effects of the vibrational leveb are implicit in the measured rate coefficient. A Boltzmann distribu- 
tion of vibrational levels prevails at the pressures where most laboratory studies of aeronomic reac- 
tions have been made. Yet, above some altitude in the mesosphere or lower thermosphere, the 
vibrational states of ozone display a non-thermal population. At these altitudes the use of published 
rate coefficients, especially those invol jig appreciable activation energies, is likely to be incorrect. 

Production and Loss Processes 

Vibrationally excited ozone is formed in the atmosphere by the reaction 

0 + ©2 + M O3 + M k^ (R7) 

where M is either O2 or N2 , and by excitation in collisions with molecules in the high energy tail of 
the thermal distribution. The products of R7 have 1.08 cV, the dissociation energy of ozone, avail- 
able for distribution among excitation and translation. The experiments of von Rosenberg and 
Trainor (1974) imply that half of the exothermicity of R7 iies in vibrational modes of the ozone 


[ 03 (i)l 

[03] 
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molecule, where on the average, the and modes combined receive 1.6 quanta of energy and 
the ^2 mode contains 3 . 7 . Table 1 lists 19 vibrational states plus their wavenumben and energies 
which can be excited with a maximum of one quantum each in the and modes and up to four 
quanta in the V2 mode. For an average initial energy of 0 .S 4 eV, states hi^er than those listed in 
Table 1 are likely excited while colliaional exdtation preferentially populates the lowest lying states. 
As will be evident, the available data do not merit a detailed treatment of all states and we therefore 
adopt the simplified method described below to estimate the atmospheric abundance of vibrationaUy 
excited ozone. 

Little information is available concerning vibrational relaxation in O3, however, the laboratory 
results of Rosen and Cool ( 1975 ) allow one to divide the process into four classes which operat: on 
different time scales. The quenching of excited O3 proceeds via a rapid exchange of energy between 
nearly resonant Vj and modes, a slower transfer from Vj and >'3 states into ^2, followed by a still 
slower quenching of 1^2 modes to the ground state. The quenching of Vj and V3 modes to 03(0,0,0) 
is much less efficient than that of i»2. The quenching rate kg ■ 2.4 x 10 "^^ cm“^s"^ appropriate 
for an O2-N2 atmosphere has been reported by Rosen and Cool ( 1973 ) and is in good agreement 


with the result of von Rosenberg and Trainor ( 1974 ). Analysis of relaxation mechanisms in ozone 
by Rosen and Cool ( 1975 ) implies that this kg value refers to the quenching: 


03(0, 1 , 0 ) + M ->03(0,0,0) + M 
Processes of the form : 

(R8) 

03(0,0, 1 ) + M -> 03(0,0,0) + M 

and 

(R 9 ) 

03(1.0.0) + M -♦ 03(0.0.0) + M 

(RIO) 

are much slower. We here adopt a rate coefficient of 0. 1 kg for R 9 and RIO, although this value is 
very uncertain. Energy tran ifers of the form 

03(0.0. 1 ) + M •+03(0. 1 . 0 ) + M 
and 

(RID 
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(R12) 


03(1,0,0) + M -03(0,1.0) + M 

proceed at an approximate rate 3 kg while the near-resonant energy exchange 

03(1.0,0) + M-* 03(0,0, 1 ) + M (R 13 ) 

has a rate coefOdent of at least 54 kq (Rosen and Cool, 1975 ). Figure 1 summarizes the above dis- 
cussion by giving a schematic representation of the coupling which exists between the three lowest 
exdted levels and the ground state. Table 2 lists the adopted rate coefticients for each of the four 
classes of exothermic reaction. Rate coeffidents of endothermic processes have the same pre- 
exponential facton listed in Table 2 for the corresponding exothermic reaction plus an additional 
attenuation of the form exp(-E/kT) where E is the energy difference between the upper and lower 
states obtained from Table 1 . 

In the altitude region where a thermal distribution exists, the lowest lying states, 03(0, 1 , 0 ), 
03(0,0, 1), and 03(1,0,0) account for essentially all of the vibrational exdtation and we treat these 
states explidtly in the calculation. The scarcity of data concerning reltxation of the hi^y excited 
states produced in R 7 makes it impossible to give a detailed treatment. We here divide all vibra- 
tionally excited states higher than 03(1,0,0) into two classes denoted by 03(t'2) and 
collectively denoted by O3*. Class 03(*'2) consists of ozone in the form 03(0, V2,0) where V2 > 1 . 
Based on the results of Rosen and Cool ( 1975 ) for 63(0, 1 , 0 ), we assume that rate coefficients for 
reactions of 03(^2) have the same pre-exponential factors as the corresponding process involving 
03(0, 1 , 0 ) as defmed in Table 2 . Similarly, 03(V| ,i'^) consists of all states, 03(Vj,v2,V3), not 
included above. The adopted rate coefficients are analogous to those of 03(1,0,0) and 03(0,0, 1 ). 
This division of O3* into two classes is based on the observed behavior of the lower states, that the 
reactive properties of the v-, excited mode differ from those of Vj and 1^3. In particular, energy 
originally distributed into t'j and 1/3 is preferentially channelled into v-, and the fmal quenching to 
the ground state occurs via a Vj niode. 

The abundance of each class of vibrationally excited ozone is computed under chemical 
equilibrium conditions in the following manner. We assume that R 7 produces 03(i>.,) and 03(1' j , *^3) 
only, consistent with the large amount of energy stored in vibrational excitation (von Rosenberg and 
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Trainer, 1974). Based on the number of individual states contained in the two classes we assume 
that one molecule in five is jn-oduced as 0^(v2) with the remainder in 03 ( 1 '! , V 3 ). Each of the five 
subsets of excited ozone is ctmnected to ail others and to the ground state by collisional deactiva* 
tion and excitation with the appropriate rate coefficient from Table 2. An additional process 
included for each state is dissociation by solar visible and ultraviolet radiation based on the cross 
sections of Ackerman (1971) and irradiance from Nicoiet (1975). In practice, the absorption cross 
section for an excited molecule is shifted slightly in energy with respect to that for the m^und state, 
however, this has a negligible effect on the dissociation rate calculation. The only significant infrared 
radiative process is the 9.6ft emission by 03 ( 0 , 0 , 1) and the subsequent excitation of 03 ( 0 , 0 , 0 ) when 
the photon Ls reabsorbed. The Einstein transition frequency for the process: 

03(0,0, 1) 03(0.0,0) -*• hy ( 9 . 6 m) (R 14 ) 

is A ■ 13.47s“* (Yamamoto, 1977'. To compute the absorption rate J, in , for: 

03(0,0,0) Iw ( 9 . 6 m) -03(0,0. 1) (R15) 

we assume a blackbody radiation field, B^(T), at frequency v and temperature T. The absorption 
cross section a(v) is assumed concentrated at the single frequency Vq so that: 

0 (v) ■ — f 8(v-i>q) 
me ^ 

where 8 is the Dirac delta function and f ■ 1 .91 S x 1(T^ is the band oscillator strength (Yamamoto, 
1977). The absorption rate is then 
J » 2wfdv a{u) B^(T) 

= 'i 

\c/\me-/ cxpdwo/^’* 

^'umerical evaluation yields values in the range 1 . 1 x 1 0 "^ - 1 . 1 x 1 0 “^ s~* between the mesopause 
and stratopause. The 14.1 m emission of O^tO. 1.0) is sufficiently weak to be negligible. Finally, 
transitions of the form: 

Ojtvj .V2.V3 ) -► Ojivj .V2.V3-I) + hi^ 
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will have Einstein coeffldents comparable to that for OjfO.0, 1). However, these merely transform 
one member of the O3* class into another and are therefore not considered here. 

Numerical Results 

In altitude regions where the atmosphere is sufficiently dense, coUisional excitation and de- 
exdtation of a given vibrational level balance to an exceUent approximation and a Boltzmann 
distribution is maintained. The processes which tend to destroy this local thermodynamic equi- 
librium are the chemical formation of highly exdted ozone in R7 and the emission and absorption 
of 9.6 m radiation. Figure 2 presents the computed production rates of O3*. Above 90 to 95 km 
chemical formation of O3* in R7 is the most efficient process. This source exceeds that due to 
coUisional excitation of the g^*ound state above 70>7S km and, above 90-95 km is greater than that 
due to further exdtation of lower lying vibrational leveb. Note that the source of O3* due to the 
sum of processes 

03(0,1,0) + M-^03* +M 
03(0,0, D + M- 03* + M 
Ojf 1,0,0) + M-*- 03* +M 
exceeds that due to: 

03(0,0,0) + M- 03* ■►M 

at aU altitudes despite the smaU relative abundance of the exdted states. This is due to the strong 
coupling of 03(0,0, 1) and 03(1,0,0) to 03(vj ,*>3) and 0(03(0, 1,0) to 03(v2) where resonant 
transfers of energy are possible. Figures 3 and 4 present the production and loss rates of 03(0, 0. 1 ) 
respectively, hormation of 03(0,0, 1 ) due to quenching of higher exdted states and by colUsional 
excitation of 03(0, 1, 0) are the major sources at aU .altitudes. However, absorption of 9.6 m radiation 
upwelling from below is of the same order of magnitude as these processes above 90-95 km. 

Figure 4 shows that 9.6 m emission exceeds collisions as a loss mechanism for 03(0,0, 1) between 
90 and 95 km. Once this occurs the assumption of a blackbody radiation field, used in computing 
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the rate of R 1 5, becomes invalid. However, since R 1 5 is never the major source of OjiO, 0, 1 ), 
the use of the Planck radiation function at all altitudes causes negligible error in comparison to the 
other uncertainties inherent in the calculations. 

The final results are given in Figure S as ratios of the number density of an excited state to that 
which wmild exist if a Boltzmaiui distribution prevailed. Below 74 km thi standard ozone profile 
adopted was the rocket model of Krueger and Minzner (1976) and above 86 km ozone was derived 
from the atomic oxygen profile of the U.S. Standard Atmosphere 1976 assuming a balance between 
R7 and photodissodation. Between 74 and 86 km a smooth curve was drawn even though there 
are indications of a secondary maximum in this region (Park and London, 1974). However, for the 
pn;sent work, which considers only ratios of txdted states, the details of [ O 3 ] in the 74 to 86 km 
region are not critical. Figure S shows that below 80 km all vibrational levels are populated accord- 
ing to a Boltzmann distribution. The highly exdted states, O^*, have the greatest deviation, due to 
the increasing influence of chemical production with altitude. Between 9S and 1 10 km the popula- 
tion of these vibrational leveb exceeds that expected for collisions alone by a factor of three or 
more. Quenching of these states tends to increase the abundance of the (0, 1,0), (0,0, 1), and 
(1,0,0) levels, however, the 9.6u emission tends to counteract this. Above 95 km the abundance 
of 03(0,0, 1 ) is less than expected from collisions alone. The efficiency of the 9 . 6 m emission as a 
loss process for exdted ozone affects all vibrational levels such that above 105 km all ratios in 
Figure 5 decrease with altitude despite the efficiency of chemical production of O 3 *. Note, however, 
that below 1 10 km the abundances of 03 ( 0 , 1,0) and O 3 * remain larger than predicted by a 
Boltzmann distribution. 


Concluding Remarks 

Laboratory studies show that each ozone m - . cule produced in the earth's atmosphere has an 
average vibrational energy near 0.5 cV. Insufficient inform-iion exists with which to perform ac- 
curate calculations of the vibrational relaxation, however, the «eini-quantitative results obtained 
here show that excited ozone deviates sigmficantlv from a Boltzmann distribution the ;opermost 
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mesosphere and lower thermosphere. Since the rate coefficients of several aeronomically 
important reactions vary with vibrational excitation, further laboratory studies of these processes 
and of vibrational relaxaticm would be valuable in defining the role of excited ozone in the 
upper atmosphere. 
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Table 1 


Vibrational States of Ozone, 03(V| , V2, V3) 
VHuch Can Be Excited With Vj , V3 < 1 and V2 < 4 


> 

> 

m 

0 

Wavenumber 

(cm“b 

Energy 

(eV) 

010 

701.42 

.087 

001 

1042.10 

.129 

100 

1103.16 

.137 

020 

1402.84 

.174 

oil 

1743.52 

.216 

110 

1804.58 

.224 

030 

2104.26 

.261 

101 

2145.25 

.266 

021 

2444.94 

.303 

120 

2506.00 

.311 

040 

2805.68 

.348 

111 

2846.67 

.353 

031 

3146.36 

.390 

130 

3207.42 

.398 

121 

3548.09 

.440 

041 

3847.78 

.477 

140 

3908.84 

.485 

131 

4249.51 

.527 

141 

4950.93 

.614 
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Table 2 


Adopted Rate Coefficients for Exothermic Energy Exchange 
Between Vibrational States of Ozone 


Process 

Type 

Description 

Exampte in 
Text 

Rate Coefficient 

I 

Quenching of i >2 modes 
to the ground state 

R8 

kQ -2.4x lO-^^cm^s-* 

II 

Quenching of Vj and Pj 
modes to the ground 
state 

R9.R10 

0.11CQ 

III 

Transfer to energy from 
I'j and ^3 into i >2 modes 

R11,R12 

3 Icq 

IV 

Transfer of enei^ from 
i»j to i »3 

R13 

54 kg 
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Figure 1 . Schematic of the coupling between the three lowest vibrationally excited levels of ozone 
and the ground state. The thickness of the lines connecting the various levels indicates 
the magnitude of the associated reaction rates. 

Figure 2 . Production rates of highly excited states of ozone by chemistry and collisional processes. 
Curves refer to the production of states with energy greater than that of 03( 1 , 0 , 0 ). 
Curve a. Chemical production in reaction R 7 . 

Curve b. Excitation of 03(0.0,0) by collision. 

Curve c. Further excitation of low*lying states by collision. 

Figure 3 . Production rates of 03(0,0,!). 

Curve a. Collisional excitation of 03(0,0,0). 

Curve b. Collisional excitation and quenching of all other excited states. 

Curve c. Absorption of 9 . 6 ft radiation. 

Figure 4 . Loss rates of 03(0,0, 1 ). 

Curve a. Quenching and collisional excitation. 

Curve b. Emission of 9 . 6 fi radiation. 

Figure S. Ratio of computed number densities to those which would prevail in a Boltzmann 
distribution. 
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Figure 1. Schematic of the coupling between the three lowest vibrationally 
excited levels of ozone and the ground state. The thickness of 










Curve b. Excitation of OjCO.O.O) by collision. 

Curve c. Further excitation of iow-lying states by collision. 
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Figure 4. Loss rates of OjfO.O.l). 

Curve a. Quenching and coHisional excitation. 
Curve b. Emission of 9.6/i radiation. 
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RATIO 

Figure 5. Ratio of computed number densities to those which would prevail 

in a Boltzmann distribution. 


